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T't shapes for streamlined, nonporous, two-dimensional, transonic 

; wind tunnels. It is based upon currently available 2-D inviscid 


transonic and boundary-layer analysis computer programs. Pre- 
dicted wa] 1 shapes are compared with experimental data obtained 
from the NASA Langley 6- by 19- Inch Transonic Tunnel where the 
slotted walls were replaced by flexible nonporous walls. Com- 
parisons are presented for the empty tunnel operating at a Mach 
number of 0.9 and for a supercritical test of an NACA 0012 air- 
foil at zero lift. Satisfactory agreement is obtained between 
the analytically and experimentally determined wall shapes. 


INTRODUCTION 

The present procedure for determining streamlined wind- 
tunnel liner shapes is being developed tc design a contoured 
nonporous liner for simulating free-air transonic flow about an 
infinite aspect-ratio yawed wing. Such a liner will be required 
for testing a large-chord, laminar flow control (LFC), swept-wing 
panel, which has a supercritical airfoil section, in the NASA Ames 
12-Foot Pressure Wind Tunnel at supercritical flow conditions. 

The nonporous liner shape is to be defined by a number of stream- 
lines which have been calculated for free-air flow about an infi- 
nite aspect-ratio yawed wing and which are further displaced 
according to a calculated boundary-layer displacement thickness 
in order to account for the blockage of the viscous layer on the 
real wall. 




This numerical procedure is applicable to a two-dimensional 
tunnel and a first check of its validity is provided by a direct 
comparison of the analytically-determined tunnel wall shape with 

one determined experimentally in a streamlined tunnel experiment. 

♦ 

Barnwell and Everhart have recently completed such an experiment 
for a symmetric airfoil at zero lift in the NASA Langley 6- by 
19-Inch Transonic Tunnel. For that test, the slotted walls of 
the tunnel were replaced by nonporous flexible plates which were 
adjusted during the experiment in order to simulate free-air flow 
over the airfoil section. Thus, the tunnel-wall displacement 
required to achieve equivalent free-air flow is determined as a 
result of the experiment. Their iterative experimental/analytical 
technique is similar to that presented by Goodyer (ref. 1). 

The purpose of this paper is to show that use of currently 
available numerical analysis procedures gives meaningful results 
which can be used in two-dimensional adaptive-wall tunnels 
operating in the transonic (supercritical) flow regime. The 
sections which follow describe the present method as it pertains 
to two-dimensional tunnel applications; very briefly outline the 
Barnwell-Everhart experiment; and, finally, compare results for 
tunnel wall shapes. The present procedure has, of course, a 
general utility which is restricted by our current ability to 
calculate the viscous flow field about arbitrary configurations. 
However, continuing rapid advances in computational .nachinery and 
methods will allow for better numerical modeling of tLe more 
complex phenomena associated with strong shocks and high lift. 

— . ' 

A detailed description of the streamlining procedure and 
results obtained from it have not yet been published by R. W. 
Barnwell and J. I. Everhart of NASA Langley Research Center. 

The present analytical results are compared with their experi- 
mentally determined wall shape; we appreciate and acknowledge 
this early release of their data. In this Paper their procedure 
and results will be identified as Barnwell-Everhart. 


2 


DESCRIPTION OF METHOD 


The present, numerical method for determining contoured-wall 
shapes which minimize wind-tunnel interference at transonic test 
conditions for large models is a direct application of several 
existing numerical techniques. First, a transonic analysis code 
is used to calculate the flow field about the effective inviscid 
model shape; that is, the geometric model plus a boundary-layer 
correction on the model. Most two-dimensional airfoil codes 
routinely used throughout the aerospace industry (refs. 2 and 3, 
for example) perform this task. In order to determine free-air 
streamline locations and pressures along them, one must use results 
which were computed throughout the flow field. When one rei aces 
the bounding free-air streamlines with a contoured tunnel wall, 
then some account must be made for the resulting wall boundary 
layer. In the present method, a two-dimensional boundary-layer 
code (ref. 4) is used to calculate a displacement correction along 
each streamline which is to be replaced by a wall. The edge 
condition for this calculation is determined from the local 
pressures calculated for the flow field about the model. For 
nonporous contoured walls at transonic test conditions, this 
viscous blockage correction must be made in order to prevent the 
streamlined tunnel from choking. This is the basic idea of the 
procedure; each particular application will differ somewhat in 
detail. Subsequent paragraphs indicate what was done in order 
to model the Barr*well-Everhart experiment. 

Thp symmetric airfoil used by Barnwell and Everhart was an 
NACA 0012 at zero incidence at the freestream Mach number 0.765. 
Under these conditions, theie is a small supersonic region termi- 
nated by a weak shockwave. When strong shocks are expected to 
occur, one must use a transonic airfoil analysis code which is 
conservative in order to predict the proper streamline locations 
(ref. 5). However, those used for the present application were 
nonconservative; thus, there is the possibility that a very small 
mass increase could occur at the weak shock and disrupt the 
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deduced in order to obtain uniform pressures through the test 
section. This first wall correction was intended to get a large 
part of the wall movement accomplished before putting the model 
in place. With the model in place, four successive wall correc- 
tions were made until it was deemed that the flow simulating the 
desired free-air condition had been achieved within the accuracies 
of the pressure measurement and wall-jack adjustment. 

COMPARISON OF RESULTS 

In this section, results obtained using the present analyti- 
cal procedure for determining the flexible tunnel-wall shape are 
compared with those obtained experimentally by Barnwell and 
Everhart. For convenience in presenting the data, all dimensions 
are normalized with respect to the airfoil chord length, c. Two 
sets of experimental data have been obtained. The first set of 
experiments was made to determine the contour of the flexible 
walls which resulted in a uniform pressure in the test section of 
the empty tunnel. The second set of experiments was made to 
define the flexible-wall contour for testing an NACA-0012 sym- 
metric airfoil at zero incidence with minimal wall interference. 

Empty Tunnel Results 

I 

For the empty tunnel test, the tunnel operating conditions 

were: 

Test section Mach number, » 0.902 

c; 

Test section Reynolds number, » 2.24x10 /cm 

Re , » 

Stagnation pressure, “ 1.60x10^ Pa 

Stagnation temperature, T^ » 302K 

In order to define t.‘:e flexible-plate contour experimentally for 
the empty tunnel, pressure measurements were made on both the 
tunnel flexible and rigid walls. Perfect gas relations were then 
used along with the area relation to determine the effective dis- 
placed tunnel dimensions. Since only the narrow, flexible, 
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4.333 chords ahead of the airfoil leading edge and the calculated 
free-air flow-field perturbations there are very small. If one 
were to neglect the wall viscous effects, these inviscid stream- 
lines would be matched to the fixed-tunnel section and form a 
planar sidewall downstream of the junction; the outermost stream- 
line would form the slightly deformed lateral wall, opposite the 
airfoil surfaces. 

Boundary-layer displacement effects are calculated along 
inviscid streamlines which are to be the tunnel walls using a 
two-dimensional, finite-difference, laminar-turbulent, boundary- 
layer code (ref. 4). The boundary layer in the tunnel upstream 
of the fixed-wall/flexible-wall junction must be considered since 
it determines the starting profile at the junction. The three- 
dimensional contraction geometry of the NASA Langley 6- by 19-Inch 
Transonic Tunnel was replaced by an equivalent two-dimensional 
configuration and boundary-layer edge conditions were determined 
from the resulting streamwise area distributions. Transitional 
and fully-developed turbulent boundary-layer calculations were 
made along the surface of the curved lateral wall as well as the 
straight sidewall for these edge conditions. Downstream of the 
contraction section, differences between the four solutions were 
found to be negligible; the f ully-turbulent solution made along 
the tunnel sidewall has been used as the starting solution at the 
junction for all viscous calculations along free-air streamlines. 

This two-dimensional procedure is not expected to be very 
realistic for the streamlines very close to the corner defined by 
the intersection of the model (which spans the tunnel) and the 
wall. In the vicinity of this intersection, three-dimensional 
effects are significant. These streamlines are also those which 
are very close to the stagnation point at the airfoil leading 
edge (trailing edge also if the airfoil has a large included 
angle there). The adverse pressure gradient at this peak causes 
the sidewall boundary layer to thicken appreciably and probably 
to separate; however, the rapid expansion around the leading edge 
quickly thins it out so that the influence is locally confined. 
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More will be said about this aspect later. In any case, it is 
felt that errors in the total blockage correction at each stream- 
wise location are not too great unless the separation or three- 
dimensional buildup represents an appreciable part of the entire 
correction . 

For a tunnel with flexible (or adaptive) walls on all sides, 
the local boundary-layer displacement correction, 6*(x,y), along 
each bounding streamline is simply applied to relieve the block- 
age. In the Barnwell-Everhart experiment (as well as most other 
two-dimensional adaptive wall tunnels), only the two lateral walls 
opposite the airfoil surfaces are flexible. (See fig. 1.) Thus, 
at each stream wise location, x, the 6*(x,y) must be numerically 
integrated around the tunnel's entire local cross-sectional perim- 
eter (i.e., across all bounding streamlines) and then applied as 
a displacement of only the two narrow flexible lateral walls. In 
the present case, the flexible-wall displacement at each stream- 
wise location is composed of the sum of two parts: (1) that due 

to the departure of the outer bounding inviscid free-air stream- 
line from a straight line - the compressible blockage, and 
(2) that due to the integrated viscous displacement correction 
for all bounding streamlines which are to be walls - the viscous 
(wall) blockage. Since the wall correction actually made by 
Barnwell and Everhart was zero at the upstream fixed-wall/flexible- 
wall junction, the analytical displacement to be compared with 
their experimental result must be a relative (differential) 
displacement with respect to that at the upstream junction. 

OUTLINE OF EXPERIMENT 

The technique used by Barnwell and Everhart is similar to 
that presented by Goodyer (ref. 1). Basically, their experi- 
mental/analytical procedure iterates on the contoured-wal 1 shape 
until consistency is obtained between the near-Iield and far-field 
solutions. The near-field (interior) solution is determined 
experimentally by the wind tunnel; that is, for a given wall 
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«?etting one observes pressures along the flexible wall (with the 
model in place) at the desired free-stream conditions. An in- 
verse far-field (exterior) solution corresponding to uniform 
parallel flow at « is obtained numerically; that is, using 
the measured flexible-wall pressures one solves for the shape 
producing it. This process is iterated until the wall shape and 
pressure changes are within the experimental measurement accu- 
racies. This procedure can be looked at as an iterative solution 
of coupled boundary-value problems where the direct inferior 
problem is solved by the tunnel while the inverse exterior problem 
is solve.'G numerically. Since the pressures along the flexible 
wall are measured in the presence of viscous effects on both the 
interior tunnel wall and the model, the tunnel blockage due to 
.hese effects has automatically been taken into account in their 
procedure. 

■*arnwell and Everhart used the NASA Langley 6- by 19- Inch 
Ti ansonic Tunnel (ref. 6) in order to verify their experimen- 
tal /analytical procedure for minimizing tunnel interference for 
testing at supercritical speeds. The tunnel was modified by 
replacing the narrow lateral slotted walls with jack-supported 
flexible plates in the region downstream of the contraction 
section. (See fig. 1.) The extent of this flexible wall section 
relative to the 15.24-cra (6-inch) airfoil chord length, c, is 
-4.833 1 x/c 5 3.333. The x-origin is the center of the model 
and is located in the test section as shown by the © on 
figure 1. The first stop in their procedure was to determine 
the wall shape required to produce uniform pressure for the 
empty tunnel (i.e., no model). This boundary-layer correction 
was required because the sidewalls of the existing tunnel are 
parallel. Tunnel runs were made at several Mach numbers; they 
found that the empty tunnel correction required at ~ 0.9 
produced adequate results at the lower Mach numbers. For this 
correction, the measured pressures at the flexible wall (and 
also the rigid wall) were plotted and wall corrections were 
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df'duced in order to obtain uniform pressures through the test 
section. This first wall correction was Intended to get a large 
part of the wall movement accomplished before putting the model 
in place. With the model in place, four successive wall correc- 
tions were made until it was deemed that the flow simulating the 
desired free-air condition had been achieved within the accuracies 
of the pressure measurement and wall-jack adjustment. 

COMPARISON OF RESULTS 

In this section, results obtained using the present analyti- 
cal procedure for determining the flexible tunnel-wall shape are 
compared with those obtained experimentally by Barnwell and 

i 

Everhart. For convenience in presenting the data, all dimensions 
are normalized with respect to the airfoil chord length, c. Two 
sets of experimental data have been obtained. The first set of 
experiments was made to determine the contour of the flexible 
wall.s which resulted in a uniform pressure in the test section of 
the empty tunnel. The second set of experiments was made to 
define the flexible-wall contour for testing an NACA-0012 sym- 
metric airfoil at zero incidence with minimal wall interference. 

Empty Tunnel Results 

For the empty tunnel test, the tunnel operating conditions 

were: 

Test section Mach number, M « 0.902 

00 

Test section Reynolds number, ^ » 2.24x10 /cm 

Stagnation pressure, ■ 1.60x10^ Pa 

Stagnation temperature, T^ » 302K 

In order to define t‘:e flexible-plate contour experimentally for 
the empty tunnel, pressure measurements were made on both the 
tunnel flexible and rigid walls. Perfect gas relations were then 
used along with the area relation to determine the effective dis- 
placed tunnel dimensions. Since only the narrow, flexible. 
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lateral walls could be adjusted, all of the displacement area 
correction was made symmetrically there. Successive pressure 
raeasurenents and displacement corrections were made until the 
desired uniform pressure distribution was obtained. 

In the present numerical analysis, flexible-wall contours 
for the empty tunnel have been determined using the same general 
procedure as with a model; the streamlines are, however, the 
existing tunnel lines. In order to facilitate the numerical 
analysis, the tunnel contraction section was represented by an 
equivalent two-dimensional configuration having a width of 
15.24 cm (6 inches) with the same area distribution as that of 
the physical configuration. Edge conditions for boundary-layer 
calcula.ions have been determined using perfect-gas relations 
and th'. tunnel area distribution. Fully-turbulent boundary-layer 
calculations made along the plane rigid wall of the equivalent 
tunnel configuration were found to adequately model the correc- 
tion. In order to be consistent with the experiment, the dis- 
placement correction is presented as a relative displacement 
taken with respect to that at the upstream junction between the 
flexible and rigid walls (i.e., at x/c = -4.833). The coordi- 
nate, y/c, defining the flexible-wall contour is measured 
laterally across the rigid wall (i.e., normal to the tunnel 
symmetry plane which is to be the model symmetry plane) with 
the origin as shown y the © in figure 1. 

Comparisons of the numerically-determined flexible-wall 
contour and that found experimentally by Barnwell and Everhart 
are shown in figure 2. The experimental displacements (denoted 
by the symbols) are those at the 11 wall-jack locations; inter- 
mediate values were not measured. Since the flexible wall is a 
thin plate, intermediate values should probably be interpolated 
using a spline fit. The maximum difference between tho analytical- 
and exper imental-displacement-thickness corrections is approxi- 
mately 18 percent of the experimental displacement at x/c = 0. 

This difference appears to be excessive; however, it corresponds 
to a difference in tunnel area ratios of only 0.44 percent. Wall 
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pressures were measured using transducers having an accuracy of 
il percent of full-scale deflection which, at these empty-iunnel 
test conditions, is equivalent to a total area-ratio difference 
of 0.36 percent. 


NACA 0012 Model Results 

The tunnel operating conditions for the airfoil test were: 

Free-stream Mach number, M^ » 0.765 

5 

Free-stream Reynolds number, ^ « 2.00x10 /cm 

Stagnation pressure, ® 1.42x10^ Pa 

Stagnation temperature, = 282K 

The test model is an NACA 0012 airfoil at zero lift with a chord 
length of 15.24 cm (6 inches); it is centered on the © mark 
shown in figure 1. 

The numerical procedure used in the present analysis of the 
airfoil experiment requires a free-air solution for the inviscid 
flow about the airfoil including viscous displacement effects on 
the model, which was obtained using the methods described iu 
references 2 and 3. At the upstream junction between the flexible 
and rigid tunnel walls (x/c = -4.833), 48 streamlines lying be- 
tween y/c “ 0 and y/c » 1.583 were extracted from the free- 
air solution for boundary-layer analysis. Figure 3 shows four 
representative free-air streamline contours for flow over the 
airfoil. The streamlines are i<’\entified by the value of Y(-«»)/c 
which gives the streamline position (relative to the plane of 
symmetry at y/c ■ 0) in the undisturbed upstream flow. Figure 3 
also shows the streamwiee locations of six lateral cuts across 
the tunnel sidewalls where representative calculated results are 
presented. 

Conventional two-dimensional boundary-layer equations cannot 
be integrated into a stagnation region with severe adverse pres- 
sure gradients, so some modification must be made along stream- 
lines near the airfoil/tunnel-sidewall Junction. The procedure 
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used in the present analysis was to replace four values of the 
pressure coefficients along each of four streamlines within 
0 < Y(- *)/c < 0.056 with the corresponding four values of 
pressure coefficients for the streamline at Y(- *)/c ■ 0.061. 

This is the streamwise region '*0.57 5 x/c £ -0.49 at the airfoil 
leading edge. This procedure reduces the level of the adverse 
pressure gradient to a value slightly belo' that for a predicted 
boundary-layer separation. Since the pressure gradient down- 
stream of the leading edge region where the modified pressure 
coefficients are used is strongly favorable, the influence upon 
the downstream boundary-layer solutions along these streamlines 
is negligible. It is noted that all of the modified streamlines 
are close to the airfoil/tunnel-sidewall Junction where the flow 
is fully three dimensional. As mentioned previously, the present 
analysis is not expected to be accurate either witii or without 
modifications in the vicinity of this junction. 

Figure 4 presents the boundary-layer thickness distributions 
across the tunnel sidewall (i.e., laterally) at tb'» six stream- 
wise locations indicated on figure 3. The displacement thickness 
distributions at x/c = -2.10 and x/c = -1.07 show the upstream 
influence of the airfoil pressure field on the boundary layer. At 
upstream distances beyond x/c ft -3.5, the displacement-thickness 
distribution with respect to y/c is essentially a constant. At 
x/c = -0.57, streamlines lying below that one at Y(- ®)/c “ 0.13 
are in an adverse pressure ^.radiant and those lying above it are 
in a favorable pressure gradient. The influence of ihe favorable 
pressure gradient is clearly shown (fig. 4(a)) in the range 
0.3 < y/c < 1.2 where the displacement thicknesses are less 
than those at an upstream location. The inflection point in the 
displacement thickness at x/c ■ -0.31 is the result of different 
flow acceleration rates along the streamlines, and the difference 
in distance along the streamlines relative to their respective up- 
sticam peak positive pressure coefficient values. It is noted 
that more drastic changes la the pressure coefficients than those 
used in the analysis (to negotiate the stagnation legion) gave 
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essentially Identl'^al results at this airfoil location and down- 
stream of it. The downstream locations x/c * 0.506 and 
x/c ■ 1.07 show the influence of the trai ling-edge compression 
and a representative profile in the wake region. At the down- 
stream location, x/c ■ 1.07, the Inflection point in the 
displacement-thickness profile results from the more rapid flow 
acceleration along the streamlines closer to the corner of the 
blunt trailing edge of the airfoil (plus boundary layer on it). 
Since the inviscid ^clucions used do not properly account for the 
wake effect, these displacement correction results should not be 
considered too realistic in this region. 

Displacement-thickness distributions along the four stream- 
lines shown in figure 7. are presented in figure 5. The behavior 
of these profiles is \ery similar to the pressure coefficient dis- 
tributions along the streamlines. In figure 6, the pressure co- 
efficient distributions along the two streamlines Yf-«)/c * 0.0456 
and *'(- ®)/c * 1.5833 are shown and this similarity is readily 
a^ jnt . 

The total displacement of the flexible tunnel walls at each 
x/c location is the sum of that due to the outer free-air stream- 
line plus tiu/it due to the viscous layer on all of the tunnel walls. 
The tunnel-wall-displacement correction is made as previously in- 
dicated. Comparison of the analytically- determined flexible wall 
shape with the Barnwell-Everhart experimental data is presented in 
figure 7. Again, tne symbols give the experimental displacement 
at the 11 wall-jack locations. The maximum difference In the data 
is at x/c ■ -0.667. The difference there in the two wall dis- 
placement c^'crections is less than 8 percent and the ratio of the 
two coordinates is 1.00?. The agreement between experimental and 
analytical wall contours is regarded as very good. 
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CONCLUDING REMARKS 


The present method for determining the contours of flexible 
nonporous wind-tunnel walls has been found adequate for the test 
cases considered. However, it is not immediately clear why the 
agreement with experimental data is better for the more complex 
case of an airfoil test than for the simple empty-tunnel test. 

One possible explanation is that the experimental data were deter- 
mined using a more complex procedure for the airfoil experiment 
than that used for the empty-tunnel experiment. 

The principal limitations of the method are due tc our pres- 
ent inability to calculate certain transonic flow phenomena. 

These are: 

(a) The displacement correction on the airfoil surface does 
not properly account for the viscous influence in the airfoil 
wake . 

(b) The three-dimensional viscous flow at the junciion of 
an airfoil and the tunnel sidewalls is not accounted for. 

(c) The procedure cannot be applied to the analysis of 
airfoils having a large region of separated flow. 

(d) An inviscid transonic code in conservation form should 
be used when strong shockwaves are to be expected. 

The principal conclusions which can be drawn from these 
limited comparisons are: 

(a) The agreement with Barnwell-Everhart experimental results 
provides some degree of verification of the present numerical 
procedure for liner design. 

(b) The detailed strearowise resolution (or control) required 
at the lateral adaptive wall is seen to be paused by the response 
of the sidewall boundary layer to the pressure field of the model. 
This is, of course, most pronounced near the model. 

(c) The blockage correction due to the sidewall viscous layer 
cannet be neglected in streamlining ^onporous tunnel walls. 
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(d) Blockage corrections alone will not necessarily account 
for the local viscous displacement effects which, at more extreme 
conditions, can modify the desired flow. That is, there will be 
local 3-D effects in the 2-D tunnel. 
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NASA LANGLEY 6- BY 19-INCH TRANSONIC TUNNEL 
WITH SLOTTED WALLS REPLACED BY FLEXIBLE WALLS 
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DENOTES CENTER OF MODEL LOCATION 
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Figure 3.“ Physical location of four streamlines and six lateral cuts across the 
tunnel sidewall where representative calculated results are given. 
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